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Recently, there have been several suggestions that weak Kerr nonlinearity can be used for gener- 
ation of macroscopic superpositions and entanglement and for linear optics quantum computation. 
However, it is not immediately clear that this approach can overcome decoherence effects. Our nu- 
merical study in this paper shows that nonlinearity of weak strength could be useful for macroscopic 
entanglement generation and quantum gate operations in the presence of decoherence. We suggest 
specific values for real experiments based on our analysis. Our discussion shows that the generation 
of macroscopic entanglement using this approach is within reach of current technology. 
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Strong nonlinear effects, if available, could be very use- 
ful for the generation of macroscopic superpositions and 
entanglement and for quantum information processing. 
For example, it is well known that macroscopic super- 
positions and entanglement could be generated in strong 
Kerr nonlinear media pj. An optical quantum computer 
could be realized if strong nonlinearity was available . 
However, it is extremely hard to obtain such strong non- 
linear effects using nonlinear media. Nonlinear effects in 
existing media are extremely weak compared with the re- 
quired level for generation of macroscopic superpositions 
and entanglement or for quantum information process- 
ing. 

Recently, it was suggested that weak Kerr nonlinearity 
can still be used for production of macroscopic superposi- 
tions and entanglement 0] and for linear optics quan- 
tum computation Nemoto and Munro showed that 
weak nonlinearity could be used to generate entangled 
states and used for linear optics quantum computation 
[H, |(|, which could also be used for the generation of 
macroscopic superposition 0, 13 • A key element of their 
scheme is a cross-Kerr interaction between a coherent 
state and a single-photon qubit. The phase of the initial 
coherent state changes by a certain amount when the 
qubit is |1), while it remains the same when the qubit 
is |0). An important point of this scheme is that the 
amplitude of the initial coherent state should become ar- 
bitrarily large to make the required nonlinear strength 
arbitrarily weak. The initial coherent state can easily 
gain an adequate amount of phase shift by |1) when a 
becomes very large. This effect could be used to gener- 
ate quantum entanglement with a very weak nonlinearity. 
The other schemes based on the singe-mode Kerr effect 
by van Enk and Jeong et al. U use implicitly the 
same principle of using a coherent state of a large am- 
plitude: all these schemes |3|, La, IM 111 use initial coherent 
states of large amplitudes in order to make the required 
separation between the component coherent states in the 
phase space with weak nonlinearity. 

However, decoherence effects during the entanglement 
generation process in nonlinear media have not been in- 



vestigated in these references 0, 0, @ in spite of the un- 
certainty of validity of such techniques under a real dissi- 
pative environment. As it was mentioned, the amplitude 
of the initial coherent state must increase at the cost of 
decreasing the required nonlinear strength (or required 
interaction time in a nonlinear medium). The require- 
ment of a large initial amplitude might reduce the co- 
herence time of the evolving state, i.e., an entangled sys- 
tem with a large initial amplitude might lose its quan- 
tum coherence more rapidly than an entangled system 
with a small initial amplitude did. It is thus unclear 
whether the decoherence effects (particularly dephasing) 
could be overcome by this approach when generating en- 
tanglement and operating quantum gates. In this paper, 
to answer this question, we investigate the effect of de- 
coherence when the idea of using weak nonlinearity [f| is 
applied to Gerry's scheme to generate a macroscopic 
superposition (so-called Schrodinger cat state) in a dissi- 
pative environment. 

The interaction Hamiltonian of cross-Kerr nonlinearity 

between modes 1 and 2 is Hk — hx a i a i a >2 a '2<> where a 
(a') represents the annihilation (creation) operator. The 
interaction between a coherent state, \a)2, and a single- 
photon qubit, e.g., \tp)i = (|0)i + |l)i)/\/2, is described 
as 

CM*Mi|a>2 = e^ t / fi i(|0) 1 + |l) 1 )|a) 2 

= -L(| ) 1 |a) 2 + |l) 1 |«e ie ) 2 ), (1) 

where |0) (|1)) is the vacuum (single-photon) state, a is 
the amplitude of the coherent state, and 6 = \t with the 
interaction time t. If 9 is n and one measures out mode 
1 on a superposed basis (|0) ± \ l))/y/2, a macroscopic su- 
perposition state, |$±) = (| a) ± | — a))/ yf~M± is created, 
where M± = 2 ± 2 exp[— 2|a| 2 ]. A macroscopic entangle- 
ment can be simply generated at a beam splitter with 
such a state. 

Using the dual rail logic, where the logical qubit basis 
is defined as \0 L ) = |1) <g> |0) and \1 L ) = |0) <8> |1), the 
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FIG. 1: A schematic of Gerry's scheme Q combined with the 
idea of using weak nonlinearity [f| . 



above process can be efficiently realized. For example, as 
shown in Fig.^ Gerry's scheme can be directly linked 
to the idea of using weak nonlinearity [BJ so that weak 
cross-Kerr nonlinearity can be used with a single photon, 
a coherent state, two photodetectors and two beam split- 
ters to generate a macroscopic superposition. If detector 
A (detector B) clicks, a macroscopic superposition state 
|$_) (|$+)) is obtained at mode b in Fig. ^ Remark- 
ably, it is clear from Fig. ^ that this approach is robust 
against inefficiency of the single photon source, loss of 
the single photon and inefficiency of the photodetectors. 
Those factors will cause the photodetectors to be silent 
and such cases can simply be discarded. Therefore, these 
will only make the deterministic property of the scheme 
non-deterministic but will not affect the quality of the 
obtained macroscopic superposition state. 

The main problem here is that a very large nonlinear 
effect, i.e. a very large 0, is required to gain a large 
separation between two coherent component states. It 
was pointed out that an optical fiber of about 3, OOOfcm is 
required for 6 = n for an optical frequency of ui ~ 5 x 10 14 
rad/sec using currently available Kerr nonlinearity |9(. In 
such a case, the state after the nonlinear interaction will 
be completely decohered because of the significant losses 
in the fiber. In order to circumvent this problem, a large 
initial amplitude a can be used with a short interaction 
time t. If a is very large, the same amount of separation 
can be obtained in the phase space even though 6 (= %t) 
is much smaller than n. 

Now we consider the decoherence effects in the Kerr 
medium. The decoherence effects can be induced by solv- 
ing the master equation [T^| 



dp 

at 
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where 7 is the energy decay rate. The formal solution 
of the master equation J5| can be written as p(t) — 
exp[( J + L)t]p(0), which leads to the solution for the 



initial element |a)(/9| 
ejq>[(j + L)t]\a)(/3\ =V(t)\a)(f3\ 



= exp[--(l 



e-^)(\a\ 2 + m+a0*]\Aa){Af3\, (3) 



where A = e~ 7 */ 2 . 

We numerically assess the decoherence effects as fol- 
lows. As we have explained, a photon loss at mode 1 
(modes a and a' in Fig. ^) will simply cause the success 
probability to be less than 1. However, energy loss of the 
coherent state part (mode b in Fig. ^| in the nonlinear 
medium should be seriously considered since it will cause 
decoherence (dephasing) of the obtained state. Note that 
the average energy loss per time increases as the initial 
energy gets larger. This will cause a more rapid destruc- 
tion of quantum coherence for a large a. The decoherence 
process (T>) will occur simultaneously with the unitary 
evolution by the Kerr effect (U) of the input state by 
the interaction Hamiltonian Hr. This process may be 
modeled as follows. One may assume that U occurs for 
a short time At, and then V occurs for another At. In 
other words, U and T> continuously take turn for a short 
time in the nonlinear medium. By taking At arbitrarily 
small, one can obtain a good approximation of this pro- 
cess for a certain time t (— A At), where N is an integer 
number. Let us set A(9 = %At = ir/N, then a larger N 
will result in a better approximation. 

We now use our model to analyze the behavior of the 
coherent state interacting with a logical qubit in the Kerr 
medium shown in Fig. ^ The first beam splitter, BSl, 
and the single photon prepares the logical qubit state, 
(|0 L ) + \1l))/V2 = (\l) a \0)a< + |0}„|l) o 0/V2. The total 
initial state after BSl but before the Kerr interaction in 
Fig. [2 is 

P(t = 0) = l(\0 L ){0 L \ + \0 L )(l L \ + \l L )(0 L \ + \l L ){l L \)®\a)(c 

Let us fist consider the evolution of the second term, 
|0l)(1l| ® |a:)(a|. After time t (= A^At) in the nonlinear 
medium, it evolves to 

{v(M)U(M)} N \0 L )(l L \®\a)(a\ = C\0 L )(l L \®\Aa)(Aae u 

(4) 

where U(At)p = U K {At)pU j K (At) and 
C = exp[-a 2 (l-e-T( i / JV )) 

N 

Y, ex P [- 7 (i/iV)]("- 1 )(l - exp[I X n(t/N)])} (5) 

71=1 

where T = \/l an d ex is assumed to be real without 
losing generality. Here we have defined C as "coher- 
ence parameter" since it determines the degree of de- 
coherence (dephasing) for the resulting macroscopic su- 
perposition state. The amplitude parameter A has been 
defined to quantify the average energy loss. The sec- 
ond beam splitter, BS2, and detectors A and B in Fig. 2] 
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FIG. 2: This figure shows that our approach gives a very 
good approximation throughout the whole range of a that 
we consider in this paper. The coherence parameter \C\ for a 
separation \a.e %xt — a\ — 2ao = 6 has been plotted for a — 300 
(solid line), a = 1000 (dashed line) and a = 10000 (double- 
dot-dashed line). 



perform a measurement onto the superposed basis states 
(|0l) ± \1l))/V2. The macroscopic superposition state 
obtained by the measurement is 



P±{t) = Af±(\Aa)(Aa\ ±C\Aa)(Aae w \ 

±C*\Aae z9 )(Aa\ + \Aae l9 )(Aae v 



(6) 



where JV± are the normalization factors. State p+(t) 
(p-(t)) is obtained when detector B (detector A) clicks 
in Fig. □ If \C\ = 1 (and A ^ 0) the state p±(t) is a 
pure superposition of coherent states while if |C| = it 
is simply a statistical mixture of two coherent states. Let 
us assume that one wishes to gain 9 = n with the ampli- 
tude of the initial coherent state is ao so that a macro- 
scopic superposition state \cto) ± | — ao) (unnormalized) 
could be obtained (without decoherence) . In this case, 
an interaction time t ~ n/x is required. However, if the 
amplitude of the initial coherent state is a (> ao), the 
required interaction time t is obtained from the equation, 
\ae lxt — a\ = 2ao, which can be derived from a simple 
geometric analysis. The required interaction time is then 
t(a,a ,x) = cos -1 [l - 2al/a 2 ]/x- 

Note that the state © can be changed into a sym- 
metric form in the phase space, i.e., Aa — > S and 
Aae tS — > —8. This can be done by applying the dis- 
placement operator, D(x) — exp(ira^ +x*a), where a and 
a' are annihilation and creation operators. The displace- 
ment operation can be performed using a strong coherent 
field (an additional local oscillator) and a biased beam 
splitter. This may be required for quantum information 
processing but it does not make any essential difference 
in our discussions. 

In our numerical calculation, we have chosen N = 10 6 , 
i.e., Ad = 7r/10 6 . It is clear from Fig. [2] that this value 
gives a very good approximation for the whole range of a 
in our study because \C\ rapidly converges as N increases. 
Our numerical results can be summarized as follows. The 
first effect is that the same amount of separation between 
the coherent states is gained in a shorter time for a larger 
initial amplitude, which we already learned from Nemoto 




FIG. 3: The decrease of coherence parameter |C| against time 
t. a = 3 (solid line), a = 30 (dashed line), a = 300 (dotted 
line). Decoherence occurs faster as the initial amplitude gets 
larger, x/l = 0.0125. 
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FIG. 4: The coherence parameter |C| against the initial am- 
plitude a for the same separation (ao = 3). The decoherence 
effect diminishes as a gets larger. See text for details. 



and Munro Q . The second effect is that decoherence oc- 
curs faster as the initial amplitude gets larger as shown 
in Fig. |3 which could be expected intuitively. Our ob- 
servation is that the first effect overcomes the second one 
so that decoherence effects diminish as the initial ampli- 
tude gets larger for the same amount of separation: Fig. 01 
shows that the coherence parameter \C\ increases as the 
initial amplitude a increases for a same separation. In 
what follows we present the detailed values obtained from 
our numerical study. One may need an optical fiber of 
about 3, 000/cm for 9 — n using currently available Kerr 
nonlinearity 9]. We choose T = 0.0125 that the ampli- 
tude will reduce as A w 0.533 for 15fcm while 9 = ir is 
obtained for 3,000fcm. This corresponds 0.364dB/km 
of signal loss, which is a typical value for commercial 
fibers used for telecommunication and easily achieved us- 
ing current technology 0, 0] . Note that signal losses in 
some pure silica core fibers are even less than 0.15dB/km 
[lll | . If the required amplitude for the obtained cat state 
is ao = 3 and the initial amplitude is also a = 3 (so that 
9 = it), the amplitude parameter is A w 2.7 x 10~ 55 , i.e., 
the resulting state will be virtually the vacuum [l|j. If 
a = 300, the amplitude parameter is A = 0.45 so that 
the "effective" amplitude, \Aa — Aa 10 1/2, of the cat state 
calculated from the separation between the two coherent 
states is « 1.35. In this case the coherence parameter 
is |C| ~ 0.047, i.e., the state is almost completely deco- 
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hered. If a = 3000, the effective amplitude is « 2.76 and 
the coherence parameter is |C| « 0.43. If a — 30000, the 
effective amplitude is w 2.97 and the coherence param- 
eter is |C| w 0.91, i.e., the resulting state will be close 
to a pure macroscopic superposition state. Therefore, in 
order to obtain |C| > 0.9 for the (effective) amplitude 
\a\ w 3 of the cat state, one needs the initial coherent 
state of a = 30000. In this case, an optical fiber of only 
about 190m will be required. 

One can simply produce macroscopic entanglement us- 
ing an additional 50:50 beam splitter on the state pro- 
duced in Fig. n 13]. The state © after this additional 
beam splitter becomes 

P E ± (t) = A4 ( (3) (ft 0\ ± C% (3) (p , 1 

±C*|/3',/3')(/3,/3| + |/?',/?')(/3',/3'l) (7) 

where /3 = Aa/y/2, 13' = Aae w /V2 and \0,j3) = \/3}\/3). 
Applying displacement operators Di{x)D2{x), this state 
can be transformed to a symmetric form as 

p% (i) =JV±{\S,5) (S, 6\ ±C'\6, S) (-6, -S\ 

±C'*\-6,-6){6,5\ + | - 5,-S){-S,-6\) (8) 

where x = -{§ + (3')/2, 8 = (/? - /?')/2 and C" = 
C exp[2A 2 a 2 i sin0]. Note that the state J7J| and the state 
|jSJ contain the same amount of entanglement because lo- 
cal unitary transformations do not increase nor decrease 
the amount of entanglement. The state |(BJl can be rep- 
resented in a 2 x 2 Hilbert space by defining the basis 
as {|$+), |$-)}. The state pf in Eq. © in the 2 x 2 
Hilbert space spanned by the new basis is then 



C'*). We have numerically calculated the negative eigen- 
value A_ of the partial transpose for pf as a measure of 
entanglement [14| . In Fig. 5, one can clearly see that 
the degree of entanglement E (= — 2A_) increases as the 
initial amplitude becomes larger. 

Our discussion clarifies that an inefficient single pho- 
ton source, two inefficient detectors, and weak nonlin- 
earity, beam splitters and a coherent state source are 
required resources for generation of macroscopic entan- 
glement, and such a method can overcome the decoher- 
cncc effect. In the same manner, weak nonlinearity can 
also be used for quantum gate operations Q in the pres- 
ence of decoherence. However, it should be noted that 
some effects in the nonlinear media such as phase noise 
may not be negligible in real experiments. It is also a 
separate problem to investigate the decoherence effect in 
nonlinear media for the other schemes 0, Q with weak 
nonlinearity. 
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FIG. 5: The measure of entanglement E of the state obtained 
by the process described in Fig. 1 against the initial amplitude 




(9) 



where K = M%{2 + C + C'*), V = -M +y /M + M-(C - 
C'*), D = M + M_(2 + C + C'*), R = M+M_(2 - C - 
C'*), W = M-y/M+M-{C - C"*), Z = AP(2 + C + 



Note added — After we submitted our paper, an up- 
date to Ref. was published, where the authors 
discussed decoherence of qubits for their computation 
scheme. 
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